Pressure of massless hot scalar theory in the boundary effective theory framework 
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We use the boundary effective theory (BET) approach to thermal field theory in order to calculate 
the pressure of a system of massless scalar fields with quartic interaction. The method naturally 
' separates the infrared physics, and is essentially non-perturbative. To lowest order, the main ingre- 

dient is the solution of the free Euler-Lagrange equation with non-trivial (time) boundary conditions. 
' We derive a resummed pressure, which is in good agreement with recent calculations found in the 

literature, following a very direct and compact procedure. 
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I. INTRODUCTION 



In finite-temperature field theory it is known that a naive implementation of perturbation theory for 
the calculation of Feynman diagrams is ill-defined in the presence of massless bosons. This is due to the 
appearance of severe infrared divergences, brought about by the vanishing bosonic Matsubara mode in thermal 
propagators. These divergences plague the entire series, which becomes essentially meaningless [if. As a result, 
one is forced to resort to resummation techniques that reorganize the perturbative series, and resum certain 
' classes of diagrams, in order to extract sensible results. There are several ways of performing resummations and 

rewriting the degrees of freedom more efficiently in terms of quasiparticles; we refer the reader to the reviews 
and to Ref. Q for a discussion and a list of specific references. All such techniques are designed to 
partially tame the infrared divergences, creating a non-zero domain of validity for weak-coupling expansions. 
Nevertheless, the zero-mode problem remains, and the region of validity of resummed perturbative treatments 
can not be indefinitely enlarged. 

In a recent paper [8], we have proposed an alternative approach to thermal field theories, denoted by 
boundary effective theory (BET) . The central idea of the method is to respect the double integral structure of 
the partition function in the functional integral formalism, 

[DM*)]p\P;<h,<t>o], ( la ) 

where 

p\P;<h,<h] = J [D<t>{r,x)] e~ s M (lb) 

#(O,a>)=0G8,as)=<M«O 

is the diagonal element of the functional density matrix and S[<f>] is the Euclidean action of the theory. This 
approach, based on the calculation of the density matrix, was already used in [6| to construct dimensionally 
reduced effective actions, and in fl to study the thermodynamics of scalar fields based on a semiclassical 
approximation. The functional density matrix formulation of quantum statistics was discussed in [7J. 

As it will be clear in the sequence of this article, the double integration approach is essentially different 
from the one where a single functional integration over periodic configurations <p(r, x) is performed. In the 
BET approach, the protagonist is the quantity <j>o{x) — the field eigenvalue in the functional Schrodinger field- 
representation: <fr\(j>o(x)) — (f>o(x)\(po(x)) . Indeed, any thermal observable can be constructed by integrating 
the appropriate functional of 4>o(x) over the fields 4>q(x) weighted by the corresponding diagonal element of the 
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density matrix. The imaginary time evolution can be viewed as an intermediate step which calculates the weights 
for the effective theory of static 0o _ne lds. That effective theory defines a quantum statistical problem which 
encodes all information about thermalization; one is led to compute correlations of the field cj>o(x), connecting 
the theory to physical quantities. 

The field 4>o(x ) has still another remarkable property: it is the zero (static) component of the dynamical 
field </)(t,x). Indeed, to each dynamical configuration <p(r, x) there corresponds a static configuration <fio(x) = 
4>(0,x) = <j)(/3,x). We say that <j>o(x) is the (time) boundary value of <f>(r,x). The difference </)(t,x) — <f>o(x) 
vanishes at r = 0, /?, and so it can be expanded in a sine-Fourier series with non-zero frequencies Cu n = mr//3, n = 
1,2,.. .. From these considerations, we conclude that the effective theory encoded in p[j3] <po, <po] also contains 
all the infrared physics, and the double integral structure of Z naturally separates the potentially divergent 
modes. 

The main result of Ref. Q was the analytic calculation of the one-loop effective action for the boundary 
field in a scalar theory. Following the standard recipe, we used the saddle-point approximation, but in two steps: 
first, we kept <fio(x) fixed, and expanded the action around the classical configuration ^ c [^o] — that has <j>o(x) 
as its boundary value, and integrated over quadratic fluctuations vanishing at r = 0,/3; then, we expanded 
the resulting expression around the saddle-point <fio = 0, and took into account quadratic fluctuations of the 
boundary field. We showed that the one-loop effective action at finite temperature has the same expression as at 
zero temperature if written in terms of the classical field <p c and if we trade free propagators at zero temperature 
by their finite temperature counterparts. Besides, we explained how to obtain a renormalized effective action 
in the case of a A^!> 4 /4! theory. 

Now we address the problem of computing the pressure of the single- well quartic theory in the (problem- 
atic) massless limit and testing the aforementioned advantages of our framework. The weak-coupling calculation 
of the pressure for the massless hot scalar theory is an enterprise of about twenty years so far [9 13], with recent 
results to order g 8 log g given in Ref. 14|. Moreover, attempts to reorganize the perturbative series have followed 
different paths, generally introducing one or more variational parameters in the pressure, to be maximized in 
the end. Among these non-perturbative methods we find screened perturbation theory (SPT), introduced in 
thermal field theory in Ref. |15| , which can also be implemented in the framework called optimized perturbation 
theory (OPT) presented in Ref. 16 1, the linear 5-expansion (LDE) [l3|, and the so-called 2PI or ^-derivable 
methods [18|. All these methods, when applied to the thermodynamics of the scalar field, are remarkably more 
stable than the weak-coupling expansion (see Refs. 0, for results using SPT, Refs. [21} for LDE, and Refs. 
ii-Ei] for 2PI). 

Due to the crucial role played by classical solutions c [</>o] i n the present approach (they are supposed 
to satisfy the Euler-Lagrange equation for arbitrary values of the coupling constant), the resulting pressure 
is essentially non-perturbative. Some technical difficulties arise, though. Decomposing the field as </>(r, x) = 
0c[0o]( r i x) + r](r,x), where rj(0,x) = r)(j3,x) — 0, introduces ultraviolet (UV) divergences in the calculation. 
In order to obtain a finite result, we use the renormalized effective action as derived in Ref. [sj]. In addition, 
spurious UV divergences appear if one naively performs the saddle-point approximation and ignores higher 
order terms in the fluctuation j](t,x) and in the fluctuations of the boundary field. A finite expression for the 
pressure is obtained by taking into account the first term in the expansion of the self-interaction. The procedure 
to avoid UV divergences can, in principle, be extended, allowing for a systematic calculation of higher order 
corrections. 

The structure of the paper is as follows: in Section UH we discuss the saddle-point approximation for 
the functional integration with fixed boundary configuration c/)q(x), and how to implement corrections to that 
approximation; besides, we write the partition function in terms of the effective action for the boundary field; 
in Section IIII1 we perform the second functional integration using, again, the saddle-point method, and a 
renormalized expression for the pressure is obtained; finally, in Section [IV| we present our conclusions. 
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II. RENORMALIZED PARTITION FUNCTION IN TERMS OF THE BOUNDARY FIELD 



Let us consider the Euclidean action, 



s[</>] = [ (d i x) 1 



'-<t> 2 + U{4>) 



(2) 



where [d x) E is a shorthand for drd 3 x. In this paper, U{<f>) — A</> 4 /4! and too = 0. For single-well potentials 
like U(4>), the unique saddle-point of the action S is the trivial vacuum <f) — 0. However, in the density matrix 
approach, the functional domain of integration is partitioned in classes where all field configurations have the 
same time boundary. The restriction of S to one of those classes (for instance, the one with boundary value 
4>o (x)) has a non-trivial saddle-point </> c [0o]j solution of 



U E ct> c {x) + m 2 <f> c (x) + U'{<j> c {x)) = , 
C (O, x) = c (/3, x) = (j> Q (x) , 



(3a) 
(3b) 



where D E — — (<9 2 + V 2 ) is the Euclidean D'Alembertian operator and x denotes (r, x). The first functional 
integration, being performed over configurations inside a certain class, will be dominated by fluctuations in the 
vicinity of the saddle-point of that class. In particular, the correspondence <j>o i-> c [0o] is 1:1 f° r single- well 
potentials. Therefore, one obtains an optimized spanning of the domain of integration in the calculation of 
Z by collecting the contributions from quadratic fluctuations around a line of saddle-point configurations, as 
suggested by Fig. [TJ One is naturally led to a two-fold saddle-point approximation. It is worth remarking that 
4> c (t, x) is, in general, a non-periodic function of r in the sense that d T (f> c (0, x) ^ d T <fi c (f3, x). One can verify it 
even in the simple case of a free theory. This is another important difference between BET and other methods. 

The explicit dependence of <p c on O is not known, except in very special cases. In we obtained the 
following recursive relation for the classical solution, 



cp c [(p ](T,x) 



d 3 x' 4> (x') d T ,G Q (T, x;t',x') 

■0 



t'=0 
t'=0 



J dr' J d 3 x' G (t, x- t', x') U'(Mt',x')) , 



where Go is a Green function of the free operator, 

(D E + m 2 ) G (x, x') = 6 (4) [x - x 1 ) 
G (r, x; 0, x') = G (r, x; /3, x') = . 

For a fixed boundary configuration, the fluctuations rj around <p c vanish at r = and /3, 

</)(t,x) = ^ c (r, x) + t](t, x) 
V (0,x)= V (f3,x) = 0. 

In terms of rj, the renormalized partition function reads 



(4) 



(5a) 
(5b) 



Zr = J [DM*)] 



[Dt](t,x)} e 



-S[4> a +-q] + C.T 



(G) 



(7) 



r;(0,a;)=77(^,a!)=0 

where C.T. are counterterms to be chosen. One can think of each 
dynamics of the fields rj takes place. 

In the vicinity of (j) c , the action is approximately quadratic, 

1 

2 



as a background field around which the 



5[0 c + 77] = + - (d 4 x) E r ] (x) D E + TOq + U" ((f>c(x)) V (x) + 0( V 3 ) 



(8) 
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FIG. 1: Pictorial description of the partition of the functional domain of integration for Z in the BET approach. The 
whole domain is indexed by the the zero-mode line. A given configuration cj>(r, x) lives on the vertical fiber over the 
boundary value <f>(0, x) = x) in the bottom. Inside each vertical line, the Euclidean action is minimized by the 
saddle-point configuration (bullet). The line of periodic configurations is also shown. 



It is convenient to introduce the Green function: 

[D E + ml + U" (0 c (x'))] G[4> c ](x; x') = <5( 4 >(x - x') (9a) 
G[<j> c ] (r, x; 0, x') = G[0 C ] (r, x; /3, x') = . (9b) 

In particular, G[0] = Go, the free propagator defined in Eq. (JSJ) . The Gaussian integration formally yields 

WoN] e- s[ ^ ] + CT (det GI0J) 1 / 2 . (10) 



We should mention that the present saddle-point approximation can be good even in a strong-coupling regime. 
In fact, the classical solution is supposed to be exact for arbitrary values of the coupling constant. Besides, one 
can systematically improve the saddle-point approximation by expanding the cubic and quartic ^-interactions, 
and calculating diagrams with lines of G[(/) c ]: 

Z R \fi\ = JpM^e-^ + ^-idetGi^}) 1 /^-^ , (11) 

where 

e~ A ^ = l-^ J(d 4 x) E G 2 [<j )c ](x,x) + ^ J(d 4 x) E G[<j )c ](x,x)+0(3 loops), (12) 

with Gi being a mass counterterm. The corresponding loop expansion is often called semiclassical series. 
In we have shown that the renormalized 1-loop effective action for the boundary field </>o is: 

PTnlfo] = S[0 C ] + \ Trlog (G~ l [(j) c ]) + \ Trlog(C[0 c ]) - C.T. , (13) 



where 



C[^} = = [drd T ,G[Mi (14) 



is the 1-loop contribution from quadratic fluctuations of the field (f>Q. 



■5 



Comparing (fTT|) and 

e -0rn[M = e -S[0 c ]+c.T. (detG [ 0c]) i/2 (detC^J)- 1 / 2 , (15) 

one can write 

Z R [p] = (detC ) 1/2 J [DM*)} e-? r *^ e" s *] , (16) 

where Co = C [0] and 

&[0 C ] = A[4> c ] - i Trlog(C[0 c ]) + l - Trlog(Co) . (17) 

In this paper, the terms involving Sj and the interacting part of S[(t> c ] will be dealt with order by order in a 
loop expansion that will be discussed in the next section. 



III. PRESSURE IN THE BET APPROACH 



Our strategy is to perform the second functional integration using again the saddle-point approximation. 
The saddle-point of r^J^o] is the trivial configuration (f>o = 0. Therefore, it is natural to expand Tr in terms of 
the following n-point functions: 

°° 1 f 

r«[0o] = ^2 ^7 / r fl V 331 '- • • > a3 «) <M^i) • • • (f>o(x n )d 3 Xi . ..d 3 x n , (18) 



where 



o(a3i) • ..6<po(x r , 



Tfl (fl) • • • j "Era) — 



(19) 



3=0 



Up to quadratic order in </>o, we have (see Ref. 

/SrflfooWr^ + i / /3rg ) (a !ll aj 2 )0o(a! 1 )^ o (aj 2 )d 3 a ; 1 d 3 a: 2 , (20) 



where = -V"tt 2 /(90/3 4 ) is essentially minus the pressure of an ideal gas of free massless bosons, V is the 
volume, and 

(2nfS( Pl +p 2 ) r *fa'ft»^) = 2 ' p il tanh^ +/3m 2 (| Pl |;/?) , (21) 

with 

2 A tanh/3A:/2 / /?fc 



W*^ = ^— + ■ (22) 



The mass counterterm chosen was 



A / A A a o 
2 7 (2tt)< 



Ci= o / 7dlA° F (g), (23) 



where Ap(q) — l/q 2 is the zero temperature (massless) free propagator in 4-dimensional Euclidean Fourier 
space. 

Substituting Eq. (|20p in (1161) with Si — 0, and performing the quadratic integration over <^o, one obtains 
the saddle-point approximation for Z, 

Z sp [p] = e~' 3r H ) (detCo) 1/2 (det/3rg ) )- 1 / 2 = e? VP » . (24) 
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In Appendix A, we show that 

, N 7T 2 1 f A d 3 k ( , coth^fc/2\ 



The second term on the r.h.s. of Eq. (|2"5)l can be identified with a series of daisy-diagrams, where the 
petal is given by ([22]) . As one can easily check, the 0(A) term in P sp is UV divergent. The 2-loop diagram 
carries the divergence. We will show that such a spurious divergence is removed when we consistently add the 
remaining 2-loop corrections to the saddle-point approximation. 

We have to be careful to identify the good propagator to represent the contraction of two <j)Q fields. We 
know that such a propagator should be calculated at the saddle-point 0o = 0. However, it is not obvious if the 
interacting mass should enter or not its definition. We claim that, in order to be consistent with the one-loop 
calculation of the effective action, we must use Co defined in (|T4l at <fr c = as the propagator. 

Notice that all quantities in this calculation depend on </>o through C . Therefore, it is convenient to 
define the contraction of two fields <j) c . In 6], it was shown that the propagator Co which is associated with that 
contraction satisfies 

A F = Go + Co , (26) 

where Ap is the usual free thermal propagator, and Go is defined in Eq. (JS|). A pragmatic argument in favor 
of using lines of Co (or Co) to build 2-loop corrections is that it solves the problem of UV divergences of the 
saddle-point approximation, reproducing the correct result for the C(A) perturbative contribution. We show 
that in Appendix B. 

Finally, the renormalized pressure is given by 
„ 7T 2 A 1 , f A d 3 k f / (3m 2 {k;f3) coth/3fc/2\ /3m 2 (fc; (3) coth/3fc/2 1 

Pa "W = m~um i ~w££oJ (^Fi log l — 2 —) — 2 —)■ 

(27) 

Figs. 2-4 plot the pressure normalized by the ideal pressure as a function of g = VA obtained using 
different methods. Fig. 2 compares the pressure from BET with that from weak-coupling calculations up to 
g a log<?, according to Ref. 14]. The weak-coupling expansion already includes resummation from order g 3 on. 
We see that the BET approach is in good agreement with the most recent results. 

Figs. 3 and 4 compare BET with screening perturbation theory (SPT) calculations at two, three and 
four loops from Ref. [i^] over different ranges. In Fig. 3, the coupling g goes from to 8. It shows a complete 
mismatch between BET and the SPT curves for, say, g > 5. However, in that range the SPT curves are not 
reliable either. In Fig. 4, we add the g s log g weak-coupling curve to the comparison. Around g — 3, SPT 
curves present quite a large oscillation as one goes from two to three, and then to four loops. We conclude that 
in the range where SPT shows convergence, the curve from BET seems to behave remarkably well. 



IV. CONCLUSIONS 



The boundary effective theory (BET) was introduced as an alternative approach to quantum statistical 
mechanics. One of its main features is providing a natural separation of the zero mode (static) sector, leading 
to the construction of its effective theory which results from integrating over all the (imaginary time) dynamical 
modes. In previous work, the one-loop effective action for the zero mode had been calculated. In the present 
article, we have shown that BET is also a powerful method to attack a crucial problem of the thermodynamics 
of bosonic fields: infrared divergences. 

The different strategies currently known to deal with such IR problems rely on some sort of resummation 
of diagrams of naive perturbation theory. The present calculation of the pressure using BET has the advantage 
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FIG. 2: Pressure normalized by Pideai as a function of the coupling constant g = VA in BET and weak-coupling (see 
[l^]) formalisms. The renormalization scale is /i = 2irT. 
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FIG. 3: Comparison of the normalized pressure obtained using BET and SPT calculations at two, three, and four loops 
reported in [201 ]. The renormalization scale is = 2nT. 



of performing a highly nontrivial resummation in the scope of a natural and systematic procedure. In fact, 
the effective theory for the zero mode generated the whole series of daisy diagrams very naturally. Besides, in 
contrast with those built from thermal field theory, the daisy diagrams built in the context of BET fit quite 
well recent results in the literature obtained using rather involved techniques of screening perturbation theory 
at 4-loops. 
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FIG. 4: Comparison of BET, g s logg weak-coupling (Ref. [1J]) and SPT (Ref. [20() normalized pressures in the range 
[0, 3]. The renormalization scale is \x — 2nT. 

A distinctive feature of BET is the role played by field configurations which are not strictly periodic in 
the imaginary time. Indeed, following the double integral scheme for the partition function, we have shown that 
there is a line of saddle-point configurations dominating the functional integral and all of them, but the trivial 
one, are non-periodic. 

We saw that the separation of the the field in its static and dynamical parts could lead to problems in 
the ultraviolet limit. This technical point has already been addressed in the calculation of the one- loop effective 
action. The solution is to perform a parallel calculation of the two functional integrations in the definition of 
the partition function. That led to a prescription for the present calculation which, we hope, can serve as a 
guide to extend the method to higher orders. 
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Appendix A 

In 0] we show that Co(k) = 2fctanh(fc/3/2). Expressing logdetCo as TrlogCo, we obtain 

f A d 3 k 

logdetCo- J log [2fctanh(fc/?/2)] . (Al) 



Analogously, from (|2lj) it follows that 

(27T) 



l°gdet/3lf= T yg, log 



2k tanh/3fc/2 + /3m 2 (fc;/3) 



(A2) 
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Therefore, 



log 



(det C ) 1/2 (det/3i4 2) )- 1/2 



A d 3 k , / „ „. coth/3fc/2 
log l + /?m 2 (fc;^) 



(2tt) s 



2A- 



Using that expression in (j24l) , one obtains ([25 



(A3) 



Appendix B 



We start by expanding the interacting terms directly in 



and 



1 - ^/ [d A x) E G 2 (x,x) + ^-J (d±x) E G (x,x) . 



Contracting the fields </> c using the proper symmetry factors, we obtain 

A f , , A v / „o / n / n\ C*i 



1 



- / (d 4 x) E (g 2 (x,x)+C 2 (x,x)) +^ / (d 4 x) E (g (x,x)+C (x,x)) 



The 2-loop contribution from Z sp [j3] is obtained expanding (see Ref. |8|) 



(detGt^]) 1 / 2 ^ (det Go)" 1 / 2 



1 - - / (d^GofoaO^aO 



(Bl) 



(B2) 



(B3) 



(B4) 



and contracting the fields <p c . Collecting the 2-loop contributions for Pbet = (log Zbet)//3U, we obtain 



Pi- 



2— loop 



8f3V 



{d A x) E (g (x, x) + C (x, x))* + -Qy I (d'x), (g (x, x) + C (x, xj) . (B5) 



Using (|2"Fj1) . we have 



P: 



2 — loop 



2/3U 



D 



(B6) 



where D is a zero-temperature infinite constant which can be set to zero. Finally, using (f2"3"|) and performing 
the remaining integration, we obtain that the 0(A) contribution to the pressure is finite and reproduces the 
perturbative result 

A 



2— loop 



1152/3 4 



(B7) 
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